Quantum magnetoelectric effect in molecular crystal Dya 
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Magnetoelectric properties of a molecular crystal formed by dysprosium triangular clusters are 
investigated. The effective spin-electric Hamiltonian is derived on the base of developed quantum 
mechanical model of the cluster spin structure. The magnetoelectric contribution to the free energy of 
the crystal is calculated. The analysis reveals several distinctive features of the magnetoelectric effect, 
which are not typical for conventional paramagnetic systems at low temperatures. The peculiarities 
are explained by the chirality of the dysprosium core of the molecules. 
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1. Studying the properties of materials with new types 
of electronic structure, such as topological insulators [l], 
or new types of charge or spin ordering, for example mag- 
netic monopoles in spin ice [2 , is a tradition of condensed 
matter physics. 

From this point of view, materials with toroidal mag- 
netic ordering is of considerable interest [3]-[5]. The gen- 
der element of such ordering is the toroidal (anapole) 
moment, which is the main term in the toroidal family of 
current (or spin) system's multipole expansion [6 . The 
simplest example of a system with toroidal moment is 
a toroidal solenoid. Toroidal moment T is a t-odd polar 
vector. 

Several antiferromagnetic crystals with a toroidal type 
of magnetic ordering have already been found ^InSj • Such 
antiferromagnets reveal interesting magneto-optical and 
magnetoelectric properties. The presence of the toroidal 
moment provides a possibility of a magnetoelectric effect 
observation. According to the [H |5] , induced polarization 
P - [H X T]. 

Since recently, magnetoelectric and the related mul- 
tiferroic materials have attracted considerable attention 
focused onto both improved fundamental understand- 
ing and novel desirable applications. Challenging and 
promising visions emerged, e.g. how to switch magnetism 
with bare electric fields and thus overcome the overheat- 
ing bottleneck in microelectronic devices p!QHT2] . 

In toroidal antiferromagnets, however, it is difficult 
to unambiguously separate effective quasi-isolated ele- 
ments with toroidal moment the interaction between 
which would create toroidal moment of the substance. 
The toroidal moment in antiferromagnets is a collective 
property that is owed to the interspin exchange interac- 



tion. 

Contrary to this, there is a molecular crystal consisting 
of Dy3 metal-organic triangular molecular clusters [13], 
which are characterized by toroidal moment [HI [15] . The 
cluster's distinctive feature is the zero magnetic moment 
of its ground state, despite the certain (clockwise or coun- 
terclockwise in the plane of the dysprosium triangle) ar- 
rangement of Dy^+ ion spins. By the analogy with single 
molecular magnets (SMMs) we will call the Dys clusters 
single molecular toroics (SMTs). Spin-electric effects in 
molecular antiferromagnets are recently reviewed in [fB] . 

The presence of the intrinsic toroidal ordering in the 
Dys system owes several peculiarities of the magneto- 
electric properties which are qualitatively different from 
those of traditional materials. The investigation of the 
peculiarities is the purpose of this work. 

First (see Sec. 2), we consider the spin structure of a 
Dy3 molecule and show in what way the toroidal mo- 
ment of the molecule is formed. The toroidal state ap- 
pears to be degenerate. But (see Sec. 3) there is a pos- 
sibility of preparing the non-generated states with the 
non-zero toroidal moment by the means of an applied 
external magnetic field. 

Sec. 4 deals with the spin-electric interactions in the 
Dys molecule. The spin-electric Hamiltonian is derived 
in the Appendix to this section and the magneto-electric 
corrections to the energy levels of Dys molecule is ob- 
tained. 

Sec. 5 presents the main results of our study. The 
magnetoelectric contribution to the free energy of the 
Dys crystal is calculated and the magnetoelectric effect 
(MEE) is described in detail. The isotherms and field de- 
pendencies of the electric polarization are obtained for 
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different ranges of magnetic field and temperature. The 
similarities and differences of the MEE in the system un- 
der study and conventional paramagnets (for example, 
rare-earth ferroborates) are discussed. 

We should note that the whole Dya crystal can be 
treated as a coherent array of the Dys SMTs due to 
weak intermolecular interaction [TTl [18]. This leads to 
the fact that macroscopic response of the crystal reveals 
the quantum properties of the very single molecule at low 
enough temperatures. Therefore, the magnetoelectric ef- 
fect in Dys molecular crystal bears traits of a macroscopic 
quantum coherent (MQC) effect. 

Finally, Sec. 6 brings forward the unique possibility 
of observing a "purely" quantum non-equilibrium linear 
magnetoelectric effect in the Dya molecular system mak- 
ing the old Pierre Curie's idea of MEE in molecules come 
true. 

2. According to the experimental data |18l [19], a 
molecular crystal formed by triangular metalorganic nan- 
oclusters Dys |T3] is monoclinic with space group C2/c. 
One primitive cell contains eight Dys molecules and is 
described by the parameters a = 32.917(2) A, 6 = 
18.0213(13) A, c = 17.3693(13) A, and p = 114.522(1)°. 
The skeleton of each Dys molecule is formed by three 
Dy^+ ions, located in the apices of an equilateral trian- 
gle. The triangles are arranged in a regular way with the 
relative sides to be parallel. Similar to the SMMs, it is 
proved experimentally [13l |18] that the Dys molecules 
are very weakly coupled, thus the crystal can be consid- 
ered as "paramagnetic", which means that the crystal 
is the ensemble of non-interacting dysprosium triangles 
with strong inter-ion exchange. Because the magnetoelec- 
tric properties of the crystal are owed to the symmetry 
of the triangle, we consider the three-ion system in a Dys 
molecule. 

The ground state of a Dy^+ ion, originating from mul- 
tiplet ^i^i5/2 split in a crystal field, is a Kramers dou- 
blet that separated from other states with energy gap 
A ~ 200 cm~^ [T7[ [18^. In these works the wave func- 
tions of the ground doublet are found to be very close to 
pure states \Mj = ±15/2) relative to the local axes of 
symmetry of each ion. Wave functions of the first excited 
doublet correspond to the states \Mj = ±13/2). The lo- 
cal easy axes (z^-axes, i = 1, 2, 3) lie on a plane of the 
dysprosium triangle at an angle of cp to the bisectors (see 
fig. [T]). The most favorable configuration corresponds to 
if = 90° [17 . 

States of the Dy^+ ion system in an external magnetic 
field under the approximation of the ground doublet are 
determined by Hamiltonian [17) 
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FIG. 1: (Upper) The structure of one of the three synthesized 
triangular Dys nanoclusters. Color scheme: blue Dy, red O, 
green CI, dark grey C, and white H. The dashed lines show 
the calculated anisotropy axes on the dysprosium fragments 
and the arrows show the ordering of local magnetization axes 
in the ground state of the complex [171 [I8j . (Lower) The spin 
structure of the Dys triangular cluster and the local easy axes 
orientation in respect of the laboratory XY Z reference frame. 



stant jzz = 10.6 K and the effective gyromagnetic factor 
gz = 20.7 [T8l. The eigenfunctions of the Hamiltonian in 
eq. ([T]) are the following vectors: 
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where S = 1/2 is the effective spin of the ground doublet. 
The best fit parameters are the effective exchange con- 



where "±" and "— " are the signs of vector J projec- 
tion of i-th ion onto local z^-axis {i = 1,2,3). Functions 
^4, . . . , <I>7 are Kramers-conjugate to the <I>s, . . . , <I>o rel- 
atively. 

The spin ordering of the Dys nanocluster is character- 
ized in terms of spin chirality [18 . It is clear that the 
spins in the states (i = 0, 1, 2, 3) and their conjugates 
(z = 4, 5, 6, 7) are inversely twisted, i.e. the states have 
the opposite chirality. The natural physical quantity as- 
sociated with spin chirality in this case is the t-odd polar 
vector of toroidal moment, which corresponds to the first 
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TABLE I: Toroidal and magnetic moments 
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term in toroidal multipole expansion of an arbitrary elec- 
tric current distribution. It reads as 

T = ^ / ((jr)r - 2r2j) d^r, 

where j(r, t) is an electric current density [6 . In the case 
of the Dys cluster, then we deal with localized magnetic 
moments, the expression can be transformed to 

i 

where is the radius-vector which connects the center 
of the triangle with i-th apex (see fig. [l]). 

The dimensionless values of the toroidal moment = 
Tz/Tq and the magnetic moment m = M//i (with 
jjL = ^gj/J^B and To = hjir) in respect of the laboratory 
OXYZ frame (see fig. m) for the states (i = 0, 1, 2, 3) 
are given in the Table ^ The values of and m for the 
conjugate states {i = 4,5,6,7) are different in sign. 
Energy levels of <I>^-states are as follows: 

3 1 

Eo= E7 = --jzz, Ei = -jzz - /im^H, z = 1, . . . , 6. 

(2) 

3. In this section we consider in detail the ground state 
of molecular nanocluster Dya as well as the possibility to 
induce the toroidal moment in it in an external magnetic 
field at different values of angle (/p, which specifies the 
orientation of the local magnetization axes in respect of 
the laboratory OXY Z frame. Let the magnetic field of 
strength H be directed along a side of the dysprosium 
triangle, i.e. H = (0;i7;0). If the field is below the cer- 
tain threshold value He = H^^p) then the ground state 
of the system is degenerate state (^o,^?)- Obviously, 
the expected value of the toroidal moment in this state 
equals to zero. On the contrary, in the field strong enough 
{H > He) the ground state is not degenerate given the 
local magnetization axes do not coincide with the tri- 
angle sides. In this case, the mean value of the toroidal 
moment differs from zero, see Table [ll} It is important to 
underline that this moment is induced by the magnetic 
field, therefore its sign can be changed by reversing the 
field direction. The field and angular dependencies of the 
Dys toroidal moment in equilibrium are shown in fig. [2] 



TABLE II: Ground state and toroidal moment of the Dys 
nanocluster in the magnetic field H = (0; if; 0) 
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Thus, depending on the direction of a sufficiently 
strong external magnetic field it is possible to produce 
the states of the system with both zero and nonzero av- 
erage toroidal moment. 

4. The toroidal moment is of high importance to reveal 
the spin-electric interactions in the dysprosium cluster 
and thus to give the description of the quantum MEE in 
the whole crystal. For this purpose, we consider electric 
field E perturbations and take into account non-Ising cor- 
rections in eq. Q. The electric field perturbation Hamil- 
tonian of a Dy" ion reads 

y- = -dE + T/^'^^ (3) 

where d = — e Xli^i dipole moment of the ion, 

e is the elementary charge, = 9 is the number of elec- 
trons in 4/-shell of the Dy^+ ion. Crystal field operator 
V^^^ in eq. ([s]) contains only odd harmonics and can be 
expressed as 

K°r'' = e^atrr*rt.(^i,<p,), (4) 

where index t is odd, Ytr{0^(p) are spherical functions, 
and atr are crystal field parameters. 

The simplest form of the crystal field operator V^.,^^ in 
eq. ([3| corresponds to the first-order term t = 1. The en- 
vironment of each Dy^+ ion is of y ^ —y symmetry, thus 
the resulting inner electric field Eq ~ 10^ V/cm effective 
for the ion by the other ions of the cluster is directed 
along the relevant bisector of the dysprosium triangle. In 
this case eq. ([3| for the k-th ion can be expressed as 

V",^" = eEo J2 {4'^ sin ^ + zf^ cos ^ . (5) 

The linear on the strength of the applied electric field 
corrections to the ion energy levels arise in the second- 
order perturbation theory with small parameter || V^H/W^, 
where ||y|| is the norm of the F-operator and W is the 
energy difference between ground states and the weight 
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center of excited ion electronic configurations (typically 
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FIG. 2: The dependencies of the toroidal moment projection 
onto the laboratory OZ axis: (upper) on the magnitude of 
the external magnetic field H for different values of angle cp; 
(lower) on the (p-angle at the strong enough magnetic field 
H — 2Hu, Hu — jzz/'^lJi- The dotted lines are the envelopes 
di sin if. 



Making use of the wave function genealogical scheme 
construction and the quantum theory of angular momen- 
tum [20 we derive the expression for magnetoelectric op- 
erator of the Dys complex. The ion electric dipole mo- 
ment is induced by the magnetic field, which mixes the 
states of the ground doublet \Mj = ±15/2) with the 
states of the higher doublet \Mj = ±13/2) (the doublets 
are separated from each other by the energy A ~ 200 ~^). 
Referring the reader to the Appendix for the details of 
calculations, we give here the final expressions for mag- 
netoelectric corrections to the energy levels of the Dya 
molecular cluster, which are bilinear on E and H: 



5EI 



3A [{HyEx - H^Ey) • sin((? + {H^E^ + HyEy) ■ cos(p] , 

^~ - ^(EyHy - EM ■ COS (3(p ± I) + {H^Ey + HyE^) ■ sin (s^ ± |) 
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+ 2A 
SE"^' 



(6) 



SE^"", A = 



3 p {rfdefgjiiB 



where rfd = 0.038 nm is the value of the radial integral 
for the Dy^+ ion [2l]|22]. The numerical value of 1^41 is 
then estimated to be 3.7 • 10~^^ cm^, therefore given by 
eq. (|6| corrections SE^^^ reach the value of SE^^^ ~ 10~^ 
cm m external fields - 1 T and £; - 10^ V/cm. The 
values SEY^^ valid for the conjugate states {i = 4, 5, 6, 7) 
are opposite in sigh with respect of SE^^^ {i = 0, 1, 2, 3). 

Generally speaking, there can exist another contribu- 
tion to the cluster polarization through variation in the 



exchange interaction due to displacement of the ions 
caused by an external electric field. This electrostric- 
tive mechanism is considered in [23^ 2A_^ by example of a 
spin triangle of antiferromagnetically coupled transition- 
metal (TM) ions. Because the exchange interaction be- 
tween rare-earth /-ions is weaker than that of TM d- 
ions by a factor of hundred, the electrostrictive double- 
ion mechanism is not effective in rare-earth clusters, 
while the above- considered single- ion mechanism plays 
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the leading role in the case of rare-earth systems. 

5. In this section we consider an equilibrium MEE in 
molecular crystal Dya in detail. In order to do this, we 
write an expression for the magnetoelectric contribution 
to the free energy of the crystal 



J' rr 



where partition function Z is 

7 



Z = ^ exp 



=0 



The expressions for energies Ei and magnetoelectric cor- 
rections 8EY^^ are given by eqs. ([2| and (|6| relatively. 
Quantity N represents the number of Dya molecules per 
volume unit, which is estimated on the basis of Dya 
molecular crystal lattice parameters (see Section 2) to 
be 8.5 • 10^^ cm-3. 

An external magnetic field induces the electric polar- 
ization in the crystal (a = x, z) 



dTrr 



dE^ 



(7) 



which depends on both magnitude and direction of the 
external magnetic field vector. 

5.1. In the limiting case of high temperatures or, equiv- 
alent ly, small fields (/imH << /cT), for an arbitrary angle 

the magnetoelectrical term of the free energy takes a 
quite simple form 

:Fme = • ^ ^^^j^kT ' + h sin 2(^), (8) 

where expressions Ii = Ex{H^ — Hy) — 2EyHxHy and 
I2 = Ey{Hl — Hy) 2ExHxHy stand for the magneto- 
electric invariants of point groups and Cs relatively. 
In the case cp — 90°, eq. ([8| contains only standard phe- 
nomenological magnetoelectric invariant Ii of a triangle 
symmetry group (1^3 ). The same invariant determines 
the thoroughly investigated magnetoelectric properties of 
rare-earth ferroborates ^1 . The invariant results in 
the typical for the paramagnets polarization dependen- 
cies on the magnitude (P ~ H^, see Fig. [s^) and the ori- 
entation (see Fig. ^) of the magnetic field strength vec- 
tor. However, we should point at the unusual for the con- 
ventional paramagnets non-monotonic polarization tem- 
perature dependence (see Fig. [s]^). The main difference 
is in the direct proportionality between the polarization 
and the cluster magnetization (not the quadratical as 
could be expected from the phenomenological analysis). 

Obtained in experiment similar curves could be a valu- 
able tool to determine the parameters of the system. 
Indeed, adjusting the direction of the magnetic field so 
that one of the polarization components vanishes (say. 



Py ~ sin(2a — 2(/9)), one can easily find the position of the 
local ion magnetization axes defined by angle (f. More- 
over, knowing the temperature at which the polariza- 
tion maximum is reached, one can evaluate the exchange 
constant j from the transcendental equation 



j/kTrr 



l + 3e 



5.2. In the case of strong magnetic fields at low tem- 
peratures the polarization behavior is different from that 
predicted by the phenomenological theory. Indeed, sup- 
posed that an external magnetic field is directed along 
the laboratory OX-axis coinciding with a bisector of the 
dysprosium triangle, the polarization vector is parallel to 
the field 



eU./kT ^ 2 cosh{V3^H/kT) + 1 ' 



Py = 0. (9) 



If the magnetic field is directed along the laboratory OY- 
axis coinciding with a side of the dysprosium triangle, 
the components of the polarization changes, so that the 
polarization vector is perpendicular to the field 

-3ANH smh{2/j.H/kT) 



eJ.z/kT ^ cosh{2iiH/kT) + 2 cosh{fiH/kT) 



, = 0. 
(10) 

Fig. [4^ shows the isotherms P^ = Pa{H)^ with a = x^y, 
for the mentioned directions of the magnetic field, which 
reveals the linearity of the MEE at the strong enough 
fields. The orientational dependencies of polarization are 
also qualitatively different from conventional ones, see 

Fig.[§D. 

Figs. [3] and [4] correspond to the most probable value 
(p = 90°. However, the actual value of angle cp is to be 
specified. This can be done by measuring the polarization 
components at the fixed magnetic field and temperature. 
Fig. |5] shows the dependence of the components for a 
range of the (p values at the strong magnetic field. 

5.3. In the weak magnetic field H < Hc^ = jzz/2/j.H = 
7.9 kOe regardless of its orientation the ground state of 
the system is the degenerate state (<l>o, ^7). Spin-electric 
interaction removes the degeneration, thus making the 
ground state of the Dya complex be a doublet with ener- 
gies [see eq. ([6| at (/p = 90°] 

SE = ±{A/To){T-[Exll]). 

To remind, T stands here for the toroidal moment, 
Tz/Tq = r^, see Sec. 2. According to eq. ([7|), if the temper- 
ature tends to zero then the induced electric polarization 
is 



P = (A7V/To)[HxT], 



(11) 



which is consistent with the expression P ~ [H x T], 
derived in |41 [5^ on the base of space-time symmetry con- 
siderations. 
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H II OX, T = 4.2 



H II OX, T= 1.9 K 



H II OY, T= 1.9 K 



H II OY, T = 4.2 K 
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FIG. 3: (Upper) The isotherms of polarization projection 
Px =bif^ for the case of a weak magnetic field directed along 
the laboratory OX and OY axes, Py = 0. (Middle) The orien- 
tational dependence of polarization projections Px ^ cos 2a 



and P, 



y 



- sin 2a for the case of a weak magnetic field H ■ 



kOe (at the temperature T = 4.2 K), making angle a with the 
laboratory OX-axis. (Lower) The temperature dependence 
of polarization projection Px in magnetic fields of Hx = 1 
kOe and Hx = 0.5 kOe for the case of high temperatures 
kT » /imH. 



FIG. 4: (Upper) The isotherms of the polarization for the 
cases of magnetic fields directed along laboratory axes OX 
(the upper curves) and OY (the lower curves). The dot- 
ted lines show the linearity of the MEE in strong enough 
magnetic fields. (Lower) The orientational dependence of the 
polarization in the strong magnetic field of if = 20 kOe, 
H = if {gx cos a-\- Gy sin a) . 



Fx (H =20 kOe, T= 1.9 K) 
Py(H=20 kOe, T = 1.9K) 



H II OY 




FIG. 5: The polarization dependence on angle (p, describing 
the position of the local ion magnetization axes. 
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"Ah II OY 


H = 20 kOe, 9 = 90° 


H II 0)^^ 






FIG. 6: The temperature dependence of the MEE constant 
absolute value for the most important orientations of the ex- 
ternal magnetic field (H = 20 kOe, = 90°). 



According to eqs. (p| and (10), there is a non-linear 



MEE in the case of weak magnetic fields but an experi- 
mental observation of the effect seems to be difficult due 
to small interlevel energies 5E'^^ ~ 10~^ cm~^. Indeed, 
the value dP/dH is less than IQ-^^ - IQ-^^ if < 5 kOe 
and T - 1 K. 

On the contrary, if the field exceeds the threshold value 
Hc2 = jzz/V^I^H = 9.1 kOe then the ground state is one 
of the nondegenerate states ^i...^^ also independently on 
the field orientation, except the fields of the levels cross- 
ings. In this case, the MEE is observable in experiment. 
As it follows from Eqs. ([9| and ([lO|, induced polarization 
P = aH, a ~ 10~^ — 10~^ dependent on the external 
magnetic field orientation, see Fig. [6) 

As for the intermediary region < H < Hc^^ the 
ground state can be both degenerate {^q^^j) and non- 
degenerate = 1...7) dependent on the field orien- 
tation. This results in the surges of the polarization, see 

Fig. [3 

6. The quantum properties of matter at low tempera- 
tures are always of special interest. It is remarkable that 
the peculiar quantum properties of dysprosium molecular 
clusters can occur in the macroscopic MEE, because they 
are very weakly coupled, which means that the whole 
crystal can be treated as an ensemble of particles with 
zero magnetic moments but nonzero toroidal moments. 

First of all, the fundamental possibility of a linear mag- 
netoelectric effect induced by an electric current in the 
fields H < He is worthy of being noted. Indeed, the 
ground state of a molecule is degenerate unless the mag- 
netic field H does not exceed the threshold value He. This 
means zero expected value of the toroidal moment and 
no linear MEE. The external current jz removes the de- 
generation due to interaction with the toroidal moment 



{Wi, 



-jT) producing a state with non-zero toroidal 



moment and thus resulting in linear MEE. 



FIG. 7: The plot of the orientational dependencies of the 
polarization vector projections Px and Py in the interme- 
diary magnetic field H = 8.1 kOe (T = 0.1 K), again 
H = if (ex cos a -\- ey sin a) . 



Not only the mentioned current effect but also quan- 
tum dynamical effects are interesting in the context of 
magnetoelectricity. The slow quantum dynamics of the 
toroidal moment in SMTs [27 gives an intriguing pos- 
sibility to observe non-equilibrium quantum MEE in a 
crystal of Ising-type rare-earth clusters, when the ground 
state is the doublet of states with energies (we imply 
ip = ^) E = ±(3/4)i,,(A/j;,)3 [57], where A is the 
gap between the energy levels. The wave functions of the 
states in the doublet are very close to = {\rz = 

+3)±|r^ = -3))/V^. In Dy3+ Kramers ions the splitting 
can be produced with an external magnetic field applied 
perpendicular to the plane of the ion triangle, A ~ QyH, 
where gy ^ 0.1 according to the measurements of the 
magnetization field dependencies at T = 1.9 K [18]. Thus, 
the splitting A is estimated to be 0.05 cm~^ in a 1-T 
magnetic field. 

Making use of crossed electric and magnetic fields (or 
just an electric current) it is possible to prepare the state 
close to the state characterized by = +3. In this case, 
there are oscillations of probability p{t) of finding the 
system in the state {r^ = —3) after time t [27 , so that 
the expected value of the toroidal moment depends on 
time {Tz{t)) = 3 — 6p(t), where p{t) = \{ip-\ipt)\^ with 
|Vt)=exp(-|m)|V+). 

According to eq. ( pT] ), there exists magnetoelectricity 
in the toroidal states P = 3{AN/To)[li x T], therefore 
the oscillations of {Tz{t)) result in the oscillations of the 
electric polarization P{t) = 3AN{l-2p{t))H . Note, that 
the polarization P is necessarily perpendicular to the field 
H. The frequency of these oscillations is = uo{A/jzz)^ 
[27] with = 48 GHz and can effectively be controlled 
for the dysprosium based system by the transversal mag- 
netic field through the value of induced splitting A. Given 
A ~ 0.05 cm~^, the period of the oscillations is estimated 
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to be 70 /is. Obtained from the experiment [T^ relaxation 
time r ~ 1 ms is large enough to observe at least a few 
oscillations of the polarization. 

7. In conclusion, the present analysis reveals that chi- 
rality of a triangular antiferromagnetic rare-earth nan- 
ocluster with Ising-like magnetic anisotropy leads to oc- 
currence of the toroidal moment that allows specific spin- 
electric interactions and quantum magnetoelectric effect. 
The resulting quantum structure of the nanocluster is 
rather rich and versatile for manipulating by external 
electric and magnetic fields or just by current. Although 
we have considered mainly the specific rare-earth molec- 
ular crystal, based on the Dys nanoclusters, many re- 
sults can be applicable to other Ising-like nanoclusters. 



e.g. Tba, H03 etc. The rare-earth metal-organic crystals 
might be the best candidates to observe the macroscopic 
quantum MEE and other peculiar spin-electric effects 
owing to the cluster toroidal ground state with long re- 
laxation time. This empowers observation of the quan- 
tum linear MEE in a crystal due to the MEE in a single 
molecule, thus making the old Pierre Curie's idea of MEE 
in molecules come true. 
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Russian Foundation for Basic Research (projects 10-02- 
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I 

Appendix 

The purpose of this section is the detailed and clear derivation of the spin-electric Hamiltonian of the Dya triangular 
complex and the calculation of the magnetoelectric corrections (|6| to the energy levels of the system. As mentioned 
in Sec. 4, the analysis is based on the perturbation theory. The second-order corrections read 

= E ((5|dE|e,){eHC'l5> + (5lC''|e.>(er |dE|5)) . (12) 

I' ,eif 

Here \g) are the dysprosium ion states of the ground configuration in a magnetic field (/ = 3 is the orbital quantum 
number of rare-earth ions and A/" = 9 is the number of electrons in 4/-shell of the Dy^+ ion), |e^/) are the states from 
the excited configuration with I' = / ± 1, and Wi' is the energy difference between |e^/) and \g) states. The 

splitting of the configuration levels is neglected. The configuration with = / — 1 = 2 {Wi-i ~ 10^ cm~^) is 

closer to the ground levels than the configuration with = / + 1 = 4. 

The interaction operator of the ion dipole moment with an external electric field in Eq. (|3| can be written in terms 
of irreducible tensor operators di^ (/i = 0, ±1) 

dE = ^(-i)^i;_^di^, 

where di^±i = T{dx ± idy)/V2, dio = d^, E±i = t{Ex ± iEy)/V2, and Eq = E^. 

For the wave functions \g) and |e^/) construction we use the genealogical scheme, in particular \g) = ^qmIJM) , 
where J and M are the quantum numbers of the ground multiplet of the ion, qm are the numerical coefficients given 
by the following equations 

\JM)= CI^^smsII^'SLMsMl), 

MlMs 

rSLMsM,) =Y. H Gl^L. ■ ^iX^lm ■ ClX^ • fs.L.M.^M., ' ("l)^"^ • V',«l,(Ci), 

where C^^^^ and stand for the Clebsch-Gordan and genealogical coefficients respectively, and Li, ^i, M/,^, 

Ms^ are the quantum numbers of the initial therm. 

The known ^20^ identity E„„. C'f;?^^^!"™, = ■ Sr, yields 

Y,{g\d^,\e,){e,ynr\9) = ^^1"^™ (sl (^i ® (r'Yt))^^ Iff}, (13) 

e^/ nm 
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where 

{g\{d^ ® (r*yt))n™l5) = E C'i'r,(<?Mi.|er)(eHr*l'tj5)- (14) 

The matrix elements {g\T}.q\ei') of irreducible tensor operator Tj^q = ^j^tkq{i) read 
{g\n,\e,) = Vn al^-Cl^^sMs-Gl^L, -Cf.M^.,;™, ■ClX,i/2, i^^i \tk,\l'm,) , (15) 

M,MlMs,Si,Lx,Ms-^ ,Mli ,mi,mi/ ,/i 

where {Imi \tkq \ I' mi') = C'/^/^^^,/.g^^7^^=-^, {l\\tkq\\l') is the reduced matrix element. 
Substituting eq. (p!5| into Eq ([m]) we obtain 



^,nJM nJM' nLML r^LM'^ (nSL \2 (^Iki 

^^LMlSMs ' ^LM'^SMs ' ^L^ML^lmi ' ^ L^M'^^lm[ ' y^S^L^j ' ^ (2/ + l)(2l' + 1) ' 

According to [W, the reduced matrix elements in Eq. ( p!6| ) can be written as (/||ri||/') = riif^/2U + IC'/^oio 



(/'||r^y,||0 = {r%'^/^^^^^^C[^%, where (r^),^ = (/|r^|r) are the radial integrals. 
Summing up of the first three factors in Eq. (16) over and m// results in (see [20]) 



E • • = (-l)V(2n+l)(2;' + l)C,'-_ I ' 5 P I > (17) 



where j ^ ^ ^ r 6j-symbols with even indexes n. 



From the equation {lmi\Ynm\lm[) = Cl^inmy ^^lORO we get Cl^lnm ^^^^ ^^"^ ^^s. ([TtI) and ^ we find 
that 



{rtYt))nm\g) = A{ll',tn) 



N 



5 )^(;r,tn)a„{5|y„«(J)|5), 



Aar,in) = -en,(rV^^^^of^|^ | ^ | V(2Z' + 1)(2; + 1), 

where an are the Stevence parameters. 

After substituting eq. (13) into Eq. (12), we arrive to the contribution of an ion to the MEE operator projected 
onto the space of the ground multiplet functions 

Hme = 25RE E(-l)'' • ■ "t- • ■ W^^ ■ CiT?, • a„ • F„™(J). (18) 

fitr nm 

The effective magnetoelectric operator of the dysprosium complex can be obtained if one sums up the contributions 
given by eq. (18) over all three dysprosium ions, 

n^,= E BZ^El^^Y^^iJ), (19) 

Oinmk 

where the index /c = 1, 2, 3 refers to Dy^+ ions, /3 = x, z, 

: 72 E • (^"-itr - C^iTitr) • A{ll',tn) ■ W^^ ■ an, 

trV 

-72 E • (^"-1*- + ^iTi*r) • A{ll',tn) ■ W^^ ■ an, 
i?^^ = 2 E ^atr ■ C^Slr ■ m', tn) ■ W^' ■ an. 

trl' 



TDX 



nm 
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In order to avoid analysis of excessively bulky expressions we consider the simplest form of the crystal field operator 
V^^^^ see eq. ([5|, and write the coordinates of the local axes explicitly 



e^^i = (sin(/9, -cos(/p,0), 
= (cos {if^^), sin ^ ^ 
e^3 = (-cos {if - f ) ,sin {ip 



e 

,0), e 
1)^0)^ 



(0,0,-1), 
- (0,0,-1), 
2/3 - (0,0,-1), 



= (cos sin (/?, 0), 
(-sin ((^+f),cos ((/:'+ I), 0) 
(sin ((/9- I) ,-cos ((/9- I) ,0) 



From eqs. (18) and ([5| after some calculations one obtains that 



EQa2C 



(20) 



where Qap{J) = \ {JaJ(3 + JjsJa — ^^a^J{J + 1)) ^re the quadrupole moment components of the Dy^+ ion /-shell, 
coefficient 0^2 = is the Stevence parameter, and C = {21 -\- 3){rfde)'^ /W. 

As it was already mentioned, the ground state of a dysprosium ion in the cluster is a doublet that is very close to 
I ± 15/2) (relative to the local symmetry axes). Excited doublet | ± 13/2) is well separated from the ground one by 
the gap of A ~ 200 cm~^. External magnetic field H entangles states | ±13/2) with states | ±15/2), thus making the 
ground ion states in the field to be as 



9± ) 



4 



4 



In this case, the only nonzero linearly depending on magnetic field matrix element of quadrupole moment is the 
component Qxz{J) 



(k) 
9± 



(fe) 



(21) 



Substituting eq. (21) for Qxz{J) in eq- (20), one obtains bilinear on E and H corrections to the energy levels Ei of 
the Dys molecular complex: 



5EY- 



(22) 



where i = 0, . . . , 7. Symbol ai^\i) has been introduced for the sign of a magnetic quantum number of k-th Dy^+ ion 
in the state (ai^^ = ±1). After summation over /c-index in eq. (22 ), we finally arrive to eq. ([6|. 
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